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Abstract: We study the new physics(NP) related to the recent discovered 125 GeV 
Higgs by employing an important subset of the standard model(SM) gauge invariant 
dimension-six operators constructed by the the SM Higgs and gauge fields. Explicitly, 
we perform a model-independent study on the production and decays of the Higgs, 
the electric dipole moments(EDM) of the neutron and the electron, and we take 
into account the anomalous magnetic dipole moments of muon and electron as well. 
We find that, even all Higgs decay channels agree with the SM predictions, the SM 
theoretical uncertainties provide a lot of room to host NP associated with the 125 
GeV boson. A linear relation is revealed in our numerical study that fizz — Hww and 
0-6 < Hzz,ww ^ 1-4 at 95%CL with or without the EDM's constraints. The neutron 
and electron EDM's severely constrain the relevant Wilson coefficients. Therefore the 
CP violating components in the h — > WW,ZZ channels are too small, ~ (9(1CT 5 ), 
to be detected at LHC However, we point out that even the parity of the 125GeV 
boson has been largely determined to be even in the h — > ZZ channel, one should pay 
special attention to the potentially large CP violation in the h — > 77 and h — > 7Z 
channels. This should be seriously checked in the future spin correlation experiments. 
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1. Introduction 

The high energy physics community has been exciting about the recent finding of a 
standard model(SM) like Higgs boson h, the final piece of SM we had long craved for, 
at the mass around 125 GeV at LHC [1,2]. So far, except the Higgs diphoton decay 
which is ~ la higher than its SM prediction, all other Higgs decay modes agree with 
the SM predictions within the experimental accuracies [3-20]. The diphoton excess 
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could be just the statistical fluctuation or due to some unknown systematic error, 
or it indicates the existence of new physics(NP) especially those connect the Higgs 
and gauge sectors. Since last year, there have been many discussion and speculation 
about the diphoton excess, see for example [21]. Whether the diphoton excess is due 
to new physics or not will be settled down by more experimental efforts in the coming 
years. Although the statistics is still weak, the current Higgs decay data has already 
set some limits to the NP associated with the Higgs boson. On the other hand, we 
should also bear in mind that the uncertainty of SM theoretical prediction for the 
Higgs productions range from ~ 2% for the vector boson fusion(VBF) to ~ 15% 
for the gluon-gluon fusion(GF) [22], the dominate Higgs production mechanism at 
LHC, see Table |2|. As for the decay branching ratios, the theoretical uncertainties 
range from ~ 3% for h — > bb to ~ 12% for h — >■ cc [22], see Table [3| So even the 
experimental precision can be improved to reach a few percent level in the future, 
still, we will not be able to conclude the total absence of NP associated with the 
~ 125 GeV boson. The interesting question will then be, " How much room for 
new physics is still allowed due to the intrinsic SM theoretical uncertainty?" The 
main purpose of this paper is aiming for a model-independent constraint on the new 
physics associated with the Higgs based on the current and future data. 

Meanwhile, the mass scale of exotic degrees of freedom(DOF) have been pushed 
to be greater than around (0.4—8) TeV in various scenarios with different assumptions 
[23,24]. Although a general analysis is lacking, it is now wildly believed that the 
cutoff scale, A, for physics beyond SM begins to show up, should be much greater 
than the electroweak symmetry breaking scale, v ~ 250 GeV. Motivated by the 
diphoton excess and the exotic DOF(s) which is(are) heavier than the electroweak 
scale, it is reasonable to assume that the new physics effects can be captured and 
described by the dimension-six (dim-6) gauge-Higgs operators. We assume that the 
NP decouple at A ^> v where the SM electroweak symmetry is not yet broken. Below 
the cutoff, there are only SM DOFs and the SM gauge symmetry are still honored 
by the low energy effective theory. The information and effects of NP are encoded 
in the form of the effective operators and their corresponding Wilson coefficients; 
see [25, 26] for the early general discussion. It is not new to study the Higgs physics 
by using the gauge-Higgs operators, see [27-31] for the previous studies along the 
same line and [32] for recent update including the Moriond 2013 data. However, 
some constraints have been over looked by those authors. The current analysis is so 
far the most comprehensive one to our best knowledge. We found those contributions 
which had been over looked indeed play a significant role to constrain NP. 

In this work we should perform a model-independent analysis based on a subset 
of SM gauge invariant dim-6 gauge-Higgs operators. The presence of these dim-6 
operators modified the gauge-Higgs couplings so as the Higgs decay and production 
at the tree-level. We should discuss the constraints on the relevant Wilson coefficients 
based on the LHC Higgs data and compare it to a benchmark case where all Higgs 
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decay channels sit at the predicted SM values and the SM theoretical uncertainties 
are used as the experimental errors. We should show that indeed the parasitic NP 
can live comfortably within the SM theoretical uncertainties without much upsetting 
the global fit. 

In addition to the tree-level processes, the coefficients of CP odd operators can 
be further constrained by estimating the loop-induced EDM at low energy. The 
apparent drawback of using effective operators is that we do not have a sensible 
prediction of the quantum effects at the loop level. The dim-6 gauge-Higgs oper- 
ators do not close under the RG running [30]. To proceed, we consider the class 
of new physics in which the fermion EDM operators can be ignored at A 1 and the 
divergent part of loop integral be removed by the counter terms. Practically, we use 
dimensional regularization to calculate the leading contribution of CP-odd operator 
to the fermion EDM in the modified minimal subtraction (MS) scheme. Similarly, 
just for a ballpark estimation of how the (g — 2) will limit the NP, without any better 
argument, we also assume the charged lepton (g — 2) is negligible at A. We calculate 
the loop induced muon and electron (g — 2) by using the gauge-Higgs operators at 
low energy and further constrain the relevant CP even Wilson coefficients. At the 
end, we will discuss a UV-complete toy model as an example for this assumption to 
work. But we should keep in mind that this assumption is not valid for the general 
case. 

The paper is organized as follows. In section 2, we will lay out the important 
subset of the dim-6 gauge-Higgs operators. The needed Feynman rules for later 
calculation will be summarized there. The modifications to the Higgs decay and 
production due to these gauge-Higgs operators at the tree-level will be discussed in 
Section 3. The 1-loop contributions to EDM and g-2 will be given in Section 4. In 
Section 5, we present the numerical results and some remarks from the global fitting. 
In section 6, we briefly discuss how a degenerate solution can be resolved by the Higgs 
pair production cross section. As an example of the UV origin of the gauge-Higgs 
operators and to illustrate in what condition our analysis is valid, a simple model 
with scalar color octet will be discussed. A brief summary will be given there as well. 
Finally, some technical details will be collected in Appendix. 

2. Effective Lagrangian 

As discussed in previous section, we assume that the NP associated with the SM 
Higgs production and decay could be largely captured and described by a set of 
SU(3)c x SU(2)l x U(l)y invariant dim-6 gauge-Higgs effective operators and their 

lr This working assumption is not totally unrealistic. In the UV theory, when all fermions are 
in their mass basis, there is no intrinsic 1-loop diagram contribution to the fermion EDM for the 
CP phases of two vertices cancel with each other. Here we further assume the relevant Wilson 
coefficients for fermion EDMs generated by higher loops at A or RG running are negligible. 
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Wilson coefficients [28,33] 



&NP = Y& * + * & i) (i=l, 2, 12, 3) . 



(2.1) 



The relevant effective gauge-Higgs operators are 



o 3 



o 



12 



o 2 




(2.2) 



where the dual field strength tensor is defined as F^ u = le^o-i* 1 , (F = B, W, G), 
and the others are in the standard notations. Note that we have absorbed the SM 
gauge couplings and the cutoff A into the definition of the operators. The operator 
O12 gives direct SU(2)-U(l)y gauge mixing and modified the oblique parameter 
S [34] at the tree level [33] 



when the two Higgs fields are replaced by their VEVs. The current S parameter 
bound, S = O.OOIqjq [35], already sets a strong limit on c\ 2 . Similarly, the O3 
operator yields an effective QCD 6*QCD-term and thus C3 is strongly constrained by 
the neutron EDM. Therefore, we will set C\ 2 = £3 = in our global analysis which 
will be discussed in the numerical section. 

There are three more dim-6 gauge-Higgs operators which give rise to the tree- 
level modification to the coupling between Higgs and gauge fields: 



However, these operators are severely constrained by electroweak precision tests; for 
example, the SM gauge boson mass matrix will be altered by these operators at tree- 
level. Since we are focusing on the Higgs physics at LHC and it is safe to ignore these 
operators, we confine our analysis with six out of the eight gauge-Higgs operators 
given in Eq. (|2.2j ). 

The relevant hV^{k\)V' v {k 2 ) Feynman rules are summarized in Table [l| where 
(k\, k 2 ) are the 4-momentum carried by the gauge bosons, and (/1, v) are the corre- 
sponding Lorentz indices. The two gauge invariant form factors are defined as 



AS 




C12 , 



(2.3) 



(D ll H)\D v H)B tlv , {p»lrfo n (D v H)W*, \H^D U H 



(2.4) 



S^ih, k 2 ) = k^k\ - h ■ k 2 g^ , P^(h, k 2 ) 



e a ^k la k 2f3 . 



(2.5) 
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Table 1: Feynman rules for the hV^(ki)V' u (k2) vertices, where (fcijfe) are entering the 
vertices. 



h'yy 






h-yZ 




■ vg% ~ 
I -JJC12 


hgg 




i ^c 3 ~ 


hZZ 




. 2vql ~ 

% -ff-a 4 


hWW 


A 2 ™?1 _ 
1 1\2~ C 2 


a ivgl ~ 

1 ~K2~ C 2 



Also we define 



«1 = Cl + C 2 - C i2 ~ Ci + C 2 , Oi = Ci + C 2 - C12 

a 2 = s 2 iy ( c 2 



ic 12 (l 



*w0 ) - «2w (c 2 - Ci^) 



«2 = S 2 vy ( C 2 — Cit^/ 



-ci 2 (i-t^; 



CI4 = C^/(c 2 + Citj'y + Ci 2 t^) ~ c^(c 2 + Cit^j/) , 



c^(c 2 + cit^ + c i2 t 2 



(2.6) 



where sin 6^ is denoted as sw for notational convenience. Similarly, we adopt the 
following abbreviations: cos^vy — > cw, sin 29 w — > s 2 Wi tan Qw —> t w etc. From 
Table [I], one can easily obtain the corresponding hhVV Feynman rules by replacing 
the VEV, v, by the Higgs field. Due to the non-Abelian nature of SU(2)l, the 
operators C 2 , (9 2j i 2 also give rise to the tree-level jWW coupling when both SM 
Higgs fields take their VEVs. The extra jWW interaction has been overlooked in 
the previous studies. However it gives nonzero contributions to the fermion EDM 
and (g — 2). Explicitly, after electroweak SSB, the relevant Lagrangian is 



5& 



7VV 



C2gl 



w: u W a ^ + c 2 gl 



4A 2 



4A 2 



W^W a ^ - c 12 g l92 



(2.7) 



The c 2 term is equivalent to a total derivative which has no effect in the local pertur- 
bation calculation. The c 2 term modifies the canonical normalization of the kinematic 
term of SU{2) gauge fields. Thus the c 2 -corresponding •yWW form factor is same as 
in the SM but the coupling is now the SM one times — c 2 gft> 2 /A 2 . 

For the triple ■j a (ki)W + ^ (k 2 )W~ ,x (k3) coupling, where all three momenta are 
entering the vertex, the corresponding Feynman rules for c"i 2 can be spelled out 



6 ~ % 2A 2 C ^ 



in 



(2.{ 
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3. Higgs Production and decay 



The SM Higgs production cross sections and its branching ratios have been calculated 
and maintained by "LHC Higgs Cross Section Working Group" [22]. For the later 
convenience, we collect the SM results in Table ^| and |3|. Note that the inherent 
theoretical uncertainties of SM prediction for both Higgs production and decays 
range from a few percents to about ~ 20%. 

Table 2: Cross sections of the SM Higgs boson at 125GeV [22], where WH and ZH stand 
for the associated production with W and Z respectively. The uncertainties are shown in 
percentages. 



tr(pb) 




LHC 7TeV (%) 


LHC 8 TeV (%) 


LHC 14 TeV (%) 


GF 
VBF 
WH 
ZH 


15.32l 14 9 
1.205+11 
0.5729+^ 
0.3158i£? 


19.52™ 
1.578^ 

0.6966+^ 

0.39431^ 


49.85i&jj 
4.180^g 
1.504+!i 

0.88301!;! 



Table 3: Branching ratios of different decay channel for SM Higgs boson at 125GeV [22]. 
The uncertainties are shown in percentages. 



10 2 x B(h -»■ ij) 


bb(%) 


cc(%) 


rr(%) 


77(%) 


ZZ(%) 


WW(%) 


gg{%) 


lZ{%) 


57.711J 


9 qi + 12 - 2 


6.32l^ 


0.228l^ 


2.641JJ 


21.5l 4 2 


8.57l£g 


0.154l|° 



Both CP-even and CP-odd operators affect the Higgs production and decay at 
tree level, and the two kinds of contributions do not mix. It is useful to define the 
following ratios to characterize different production and decay channels of the SM 
Higgs 

rjj -> ij) a(XY -> h) 

a v r SM (/i^y)' 1XY a^XY^h) 1 1 ' 

where and (X, Y) stand for the final states and the initial particles respectively. 
At tree-level, the gauge-Higgs operators have no effects on the SM Yukawa couplings 
so that aff = 1. Since the gluon-gluon fusion(GF) is the dominate contribution to 
the Higgs production at LHC, for simplicity we will ignore the change to the other 
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production channels due to the presence of the gauge-Higgs operators. Next, we 
discuss how the gauge-Higgs operators modify 7 S9 and the relevant cey. 

3.1 Gluon fusion 

The effective operators which describe short distance features will not change gluon 
PDF and kinematics, so we just need to focus on partonic cross section, which is 



<W#3 ->■ h) = a^m 2 h 5(s - m 2 h ) 



7T 



8m, 



-T(h — >■ gg)S(s — m 2 h ) . 



Then the ratio of cross section is 

a(gg^h) Y(h -> gg) 



7gg 



vsM(gg^h) T SM (h^gg)' 



(3.2) 



(3.3) 



In SM, Higgs decays to two gluons via heavy quark loop, and the decay width 
is well-known [36] 

G F a 2 m 3 h 



T(h gg) 



where r t = and Ai is given by 



64v^7T 3 



(3.4) 



A l2 (r)=2[r+(r-l)f(r)]r- 2 , 
arcsin 2 Jt 

f{r) = { 



l n i±yT^_ 



T < 1 
T > 1 



(3.5) 



Incorporate the CP-even and CP-odd effective operator contributions, the total 
decay width can be calculated to be 



r(fc gg) 



G F A 2 ° 3 



2 






+ 


G F A 2 ° 3 





(3.6) 



which agrees with [29]. But as had already stated, the C3 is severely constrained by 
neutron EDM, only the CP-even operator O3 is relevant to the ratio of gluon- gluon 
fusion cross section 



I99 



16V2n 2 
+ G F A 2 A,(T t ) C3 



(3.7) 



3.2 Diphoton decay 

The h — > 27 decay width relates to the amplitude 

1 



327T rrih 



\M^\ 2 , 



(3.8) 
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where 



M" = + 2 -^L ai S^(k u k 2 )e\ a el, + 2 ^~a 1 P^(k 1 , k 2 )e* la e* 2 , , (3.9) 

with the SM contribution [36,37] 



(3.10) 



where Ai is from top quark loop as given in the GF section, and Ai is from the W 
boson loop, 



A^t) = - [2r 2 + 3r + 3(2r - l)/(r)] r~ 2 , 



(3.11) 



and parameter Tj = ^fc. It is straightforward to calculate the ratio of decay rates 



a 



77 



\M(h-> 71 )\ 2 
\M S M(h ->■ 77) | J 



(3.12) 



8v^7r 2 ai 



G F A 2 [^i(r w ) + |Ai(r t )] 



7r ai 



G' F A 2 [A 1 (r^) + |Ai(r t )] 



3.3 h ->■ 7Z 

The decay rate for ft, — > 7Z in the rest frame of Higgs boson is 



r = -^—(i-^\ l-^ 7Z l 2 - 



(3.13) 



The amplitude including high dimensional operators' contributions is 



JW Z = Mf^ Z) + V j^a 2 S^e^e* 2u + ^a 2 P^ Q e^ , (3.14) 



A 2 



and SM contribution [36, 38, 39] is given as 



M^ Z) =M 2 S^6l u e* 2l/ , 



M 2 = ^-(V2G F y (A F + A W ), 

A F = — <2 ^—A 1 f{r u A t ) , A w = — —A1 z (t w , \ w ) , 

C\V S W 2 S W 



(3.15) 
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where Tj = A» = (i=t, W), and fermion vector coupling Vf = 27| — 4Q/s^. 
The relevant functions are 

Af(r,A) = [/ 1 (r,A)-/ 2 (r,A)], 

2 

(r, A) = c w {4 (3 - t 2 w ) J 2 (r, A) + [(1 + 2r)t 2 w - (5 + 2r)] Ji(r, A)} , 



/i(r,A) 
/ 2 (r,A) 



+ 



2(A-r) 2(A-r) 2 

[/(T) - /(A)] 



[/(r) " /(A)] + 



A 



(A-r)' 



[<?(r) - <?(A)] , 



v^ 1 

Vl-T- 



1 arcsin ^/r 



In 



i+\/T- 



l-VT- 



l-K 



T < 1 
T > 1 



(3.16) 



Then the ratio of decay rates is 
\M{h^ 1 Z)\ 2 



a 



\M SM (h^>yZ)\< 



G F A 2 {A F + A W ) 



8y/2- 2 * 



vr a 2 



G F A 2 (A F + A w ) 



(3.17) 



3.4 h -> WW*,ZZ* 



Unlike the cases of h — > gg, 77, 7Z which are 1-loop processes in SM, the leading 
contributions to h — > WW* , ZZ* start at tree-level in SM. Moreover, the SM tree- 
level hV^V v vertex is ~ g^ v which mixes nontrivially with the form factor when 
squaring the amplitude. So aww,zz can not be expressed as |1 + ec even \ 2 + \e'c dd\ 2 , 
where e, e' are some small numbers, as in the previous cases. To obtain the h — > 
VV*(V—W,Z) decay, one needs to perform a straightforward tree-level calculation 
and take care of the phase space integration of the 3-body final state. When the 
decay final states of gauge boson V* are massless, the expression can be largely 
simplified but still not very illustrating. The details can be found in the appendix. 
With the presence of the gauge-Higgs operators, the ratio of h — > VV* decay rates 
are 



otww — 



Rx{m w jm h , 



otzz = 



Ri(m w /m h ,0) 
R 1 (m z /m h ,a z 



„ 2 32mlm 2 v R 2 {m w /m h ) 
+ c 2 



A 4 



~r 04" 



Ri{m w /m h) $) ' 
R 2 (m z /m h ) 



R 1 (m z /m h ,0) 4 A 4 ^KK.O) 
where the relevant CP even NP is encoded by two parameters: 



(3.18) 



aw 



8m w 
A 2 



-c 2 , a z 



8m|, 
A 2 



0,4 



Our results for a ggninZ agree with [29] and the explicit expression for aww,zz 
are however new. Now we are fully equipped for the later numerical analysis on the 
LHC Higgs decay data. 
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3.5 Numerical expressions for a's 

Taking A = 1 TeV as a reference point, numerically, we have 

a gg = Igg = (! + 13.90c 3 ) 2 , 
a 77 = [1 - 1.48(0! + c 2 )] 2 + 2.18(ci + c 2 - 5 12 ) 2 , 
a yZ = [1 + 0.46C! - 1.54c 2 ] 2 + 0.71 [0.555! - 1.82c 2 + 0.64g 12 ] 2 , 
a ww = 1 - 0.019c 2 + 1.3 x lO^c 2 . + 5.2 x 10~ 5 c 2 2 , (3.19) 
a zz = 1- 0.86 x 10" 3 C! - 0.0096c 2 + 2.6 x 10~ 7 c 2 + 5.8 x 10" 6 cic 2 + 3.0 x 10~ 5 c 2 
+1.3 x 10~ 5 (c 2 + 0.09ci + 0.28ci 2 ) 2 . 

By re-scaling q — > C{ x (A/lTeV) 2 , one can easily obtain the numerical expression 
for the cutoff scale which differs from ITeV. Since the h — > 77, Zj, gg are the loop 
precesses in SM, these modes are more sensitive to the NP. One can see that the 
prefactors associated with c's and c's in a ggrnri z are all around a fews. On the other 
hand, the leading contributions for h — > WW, ZZ begin at tree-level in SM. The 
NP modification to these modes are relatively two orders weaker, by a loop factor 
~ 1/167T 2 , compared to h — > gg, , y-f,'-fZ. 

4. EDM and g-2 




Figure 1: The 1-loop contributions to electron and quark EDM's and charged lepton g-2. 
The gray bulbs represent the effective gauge-Higgs operators. Note that the mirror image 
diagrams are not displayed. 

First, we set our convention. The electric dipole moment (EDM) of a fermion, 
df, is defined by the low energy effective Lagrangian 

=%DM = ~dfU(p 2 )h5C r fiuU(p 1 )F^ , (4.1) 

where u and u are the spinor wavefunction. It is well known that the fermion EDM 
starts at least at 3-loop level for quarks and 4-loop level for leptons in SM. However, 
with the CP-odd operators in the effective Lagrangian (|2.1| ), the fermion EDM can be 
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generated at one-loop level, and there are three possible contributions from the CP- 
odd hjj, h^Z, and jWW interactions, see Fig.|l| We note in passing that quarks can 
acquire nonvanishing chromoEDM if £3 7^ 0. Similar consideration of using the EDM 
to constrain the CP-odd Wilson coefficients can be found in [27,28], where only the 
contributions from Fig.|l](a) has been taken into account. As discussed in the intro- 
duction, we perform the 1-loop calculation by dimensional regularization and extract 
the finite part by MS scheme. Our result for the complete 1-loop contributions is 
summarized as 



dt 



arrif 



7T IT 



Q f aiKi(K,m h ) + a 2 2 f n2 ^ f — K 2 (A,m z ,m h ) + ^-Ki(A,m w ) 



S W S2W 



Ml' 



(4.2) 

where the definition of ai, a 2 have been given in Eq. (|2.6|) , If is fermion's isospin, and 
in our convention, e = |e| and Q e = — 1. In the bracket, the first, second, and third 
term represent the contributions from Fig.|l|(a),(b), and (c) respectively. In the limit 
that rrif <C m/j, the functions K\p take the form 



Ki{fi,x) = — 



K 2 {n,x,y) = — 



3 1 ^ 

4 + 2 ln ^ 



1 (x 2 In ^ — V 2 In ■ 



(4.3) 



where /x is the dimensional regularization scale. Similar result for the /177 contribu- 
tion can be found in [28], our result is different from [27] up to a finite term. For 
neutron, we adopt the QCD sum rule estimation 



d n = {\± 0.5) l.le(d d + 0.5d u ) + lA{d d - 0.25d u ) 



(4.4) 



given in [40] to relate the quark EDM and neutron EDM. We take m u = 2.3MeV, 
rrid = 4.8MeV, and a = 1/128.0 for our numerical study. The quark chromoEDM's 
are ignored due to the 6qcd constraint. In terms of the Wilson coefficients, we obatin 



d e = (7.00ci + 7.39c 2 - 16.07c 12 ) x 10~ 2b ecm, 
d n = (1.91ci + 10.045 2 - 16.25c 12 ) x KT 25 ecm. 



(4.5) 



From the latest bounds: \d e \ < 1.05 x 10~ 27 e cm(90% C.L.) [41], and \d n \ < 2.9 x 
10~ 26 e cm (90% C.L.) [42], we obtain two inequalities 



|66.69£i + 70.37£ 2 - 153.015i 2 | < 1 (1.6a) 
(1 ± 0.5) x |6.60ci + 34.645 2 - 56.025 12 | < 1 



[1.6a) 



(4.6) 



When combining the above two conditions together, the allowed region is a solid 
tube with tiny parallelogram cross section passing through the origin along the 
{0.735, 1.478, 1.0} direction in the {£1, c 2 , 5i 2 } space. 
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The same 1-loop diagrams in Fig. ^] with CP-even operators give rise to a/, the 
anomalous magnetic dipole moment (AMDM) of charged leptons. The AMDM can 
be extracted from the charged lepton's photon form factor, 



iM. = ieu(p 2 



u(pi), a f = y -—=F 2 (0), (4.7) 



where the momentum transfer q = p 2 — p\. Again, we perform the 1-loop calculation, 
in the unitary gauge, with dimensional regularization and in the MS scheme. We 
obtain the following result 



arrif 
a f = 2 

7T V 2 



2a 1 Q f K 1 (A,m h )+2a 2 * If 2 Q/ ^ i( 2 (A,m,,m t ) + 5 



' (4-8) 

The first and second term in the bracket are the finite part of Fig.^(a) and (b) 
respectively. The third term, which has no divergence, stems from the effective dim- 
4 (WW^) interaction introduced by operator 2 . The contribution to a/ from this 
new interaction can be obtained by multiplying the factor — [c 2 g\v 2 /A 2 ) 

to the well known SM WW'j contribution [43]. Numerically, we obtain Aa M = 
1.61(d + 0.41c 2 ) x 1(T 10 for A = 1 TeV. 

For charged lepton anomalous magnetic diploe moment, the deviation of exper- 
imental measurement from the SM prediction are 

Aa„ = a e * p - af 1 = (2.39 ± 0.79) x 1(T 9 (la) [35,44] , 

Aa e = a e * p - a s e M = -10.6(8.1) x 10~ 13 (la) [45] . (4.9) 

Since the gauge-Higgs operators are flavor blind, the resulting Acif scales as rnj so 
Aa e = (m e /m /i ) 2 Aa /J . Assuming that Aa e ^ are solely attributed to the gauge-Higgs 
operators and using the latest data, Eq.([4.9|), we obtain a best fit at Aa^ = 2.37x 10~ 9 
from the least square fit, and 

C! + 0.41c 2 = 15.25 ±5.09. (4.10) 

The allowed region is basically a wide infinite strip away from the origin on the c\-c 2 
plane. 



5. Numerical analysis 

To quantify how each Higgs decay channel differs from its SM prediction at LHC, 
we follow [46] and use the signal strength parameter fx which is defined as 

A a{X -> h)B(h -> if) 
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where X stands for the initial partons in proton, and i, j represent the decay products. 
The branching fractions are given by Bij = B(h — > ij) = Tij/T tot , and r to t is 
the actual total Higgs decay width. The r tot is related to the SM prediction by 
Ttot = Ctotrfot 1 - Then, in terms of a's and B's, C to t can be expressed as the sum of 
all contributions from the major decay channels 

t^tot ~ &bb + B tt + "cc + "77^77 + a gg&gg + a WW&WW + a ZZ&ZZ + a ZjD Zry 

~ 0.67 + 0.0023a 77 + 0.0867 M + 0.026a zz + 0.215a ww + 0.0015a Z7 , (5.2) 

where the SM branching ratios for a 125 GeV Higgs are adopted from [22] and note 
that ajf = 1 at tree-level since the gauge-Higgs operators do not modify the Yukawa 
interactions. Due to the parton distribution function, the loop induced gluon fusion 
is the dominate production mechanism (~ 87 — 88% ) at LHC for a Higgs with mass 
around 125 GeV. And among all production channels, it is most sensitive to new 
physics. Therefore, it is a fairly good approximation to take a(gg — > h)/a(gg — > 
h)sM ~ Igg- Therefore we have 

ho ~ • (5-3) 

As the common practice, a function x 2 is defined for multi-parameters fitting 

X 2 = £ iA ^ 1 . (5-4) 

i ' 

where /ij is the mean experimental value, <5j is the combined uncertainty from both 
the experimental and theoretical sides, and fii is the theoretic prediction either from 
a specific model or determined by the effective gauge-Higgs operators. We use the 
most updated Higgs data collected and analyzed by ATLAS and CMS groups, see 
Tab.[|. Since we explicitly take a(gg —> h)/a(gg —> h)su ~ Igg — a gg, i n the CMS 
h — > ZZ* channels, we only use the gluon-gluon fusion data. We drop from 
the global fitting for both groups have large uncertainties and both are via the VH 
production. 

5.1 Pseudo Global fit of ATLAS and CMS results 



To proceed, we first treat the five a's in Eq.( |5.2j ) as mutually independent free pa- 
rameters and conduct a pseudo global fit by using the current LHC data listed in 
Tab.§. The best fit locates at 

7 Jf C = 0.99, = 1.26, 0^ = 0.84, a^f = 1.30, a^ c = 1.00 (5.5) 

with minimum Xo — 4.94, and the la boundary corresponds to the contour of Xi a = 
Xq + 5.89 for a five parameters fit. And the la model-independent correlations 
between and 7 99 are displayed in Fig Basically the message from doing this 
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Table 4: The signal strength [x from Moriond 2013. 





CMS 


ATLAS 


zz* 

WW* 

77 

lZ 


0.91±gJ2 (inclusive) 
1.0i|| (qqH,VH) 
0.9i°j (GF) [13,17] 

0. 761g|J (inclusive) [13, 18] 

0.78l°;i (MVA) 

1. niolo ( cut based ) [ 14 ,i6] 

< 9 [20] 


1.7±o 5 4 (inclusive) [3,9] 

1.01 ±0.31 (inclusive) [4,10] 
1.65±{j:3o (inclusive) [3,8] 

< 18.2 [3,7] 


bb 

TT 


HCP12 : 1.31^ (VH) [6] 
l.ll^ (inclusive) [15,19] 


1.09 ±0.20 ±0.22 (VH) [5,12] 
0.7 ±0.7 (inclusive) [6,11] 
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Figure 2: The la correlations between: (a) a 77 and j gg , (b) oty/w and ~y gg , (c)azz and 
7 99 , (d)a 7 z and j gg . The best fit locations are marked by dots. 



trivial exercise is clear and simple: (1) the current LHC data more or less agrees 
with the SM, and (2) the larger the Higgs gluon-gluon fusion production the smaller 
the Higgs decay widthes and vice versa. We see that the global fit prefers a SM-like 
Higgs gluon fusion production. However, this still introduces a degenerate solution 
to C3, (1 ± 13.92C3) = ±y/1.0 ± 0.7. And the two corresponding la allowed region 
for c 3 are e [-0.167,-0.114] and e [-0.029,0.023]. This degeneracy for c 3 also 
inevitably exists in the later global fitting. We will discuss how to lift this degeneracy 
phenomeno logically in Sec. 6. 

5.2 Limits on the gauge-Higgs Wilson coefficient by using the Higgs data 
alone 

Next, we do the global fit by using the current LHC data in terms of 01,2,3 and 61,2,12- 
For a 6-dimensional parameter fitting, the 95%CL (68%CL) contour corresponds to 
^95(68) = Xmin + 12.59(7.01). To have an idea how the Wilson coefficients will be 
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constrained when the experimental sensitivity are improved and eventually compat- 
ible with the SM theoretical uncertainties, we also consider a 'fake' data set that all 
signal strengthes equal one, ft = 1.0. And we take the combined SM theoretical un- 
certainties, <jsm = \j &br + a GF-> fr° m Tab.0 and Tab.[|, as the 'experimental' errors. 
The 95%CL results are shown in Fig|| 




Figure 3: The 95%CL allowed region in the Wilson coefficient space by using the LHC 
data (blue/darker) and the SM predictions (light brown/lighter). The best fit location 
is shown at the dot(cross) for LHC data (SM). The minimum \ 2 = 6.413(0.0) for LHC 
data(SM). In subdiagram(b), we only display the SM like allowed region for C3. 

Some features of our results: (1) The CP even Wilson coefficients are well con- 
strained by the current LHC data. The LHC limits on Ci,C2 and C3 are basically 
compatible, around two times bigger, with the benchmark experimental sensitivity 
by using the SM theoretical predictions and uncertainties as input. (2) On the other 
hand, the constraints on the CP-odd Wilson coefficients £1,2,12 are poor by using 
the Higgs decays data alone. However, they seem to fall on a long line segment in 
the parameter space. This linear relation can be easily understood as roughly the 
solution to make the CP odd contributions vanish simultaneously in a 77 and a 1 z- I 
Eq.( p,19| ), since the two are most sensitive to the presence of CP-odd contributions. 

The most important message from this drill is: there still is the allowance for 
NP with Ci,2 ~ a few and c's ~ (9(100) for A = ITeV hiding in the SM theoretical 
uncertainties. 

Interestingly, from our numerical study, we find a linear correlation among oiww 
vs ctzz and the signal strength fiww vs f^zz, see Fig.f|. Our results show that at 
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95% CL azz ~ l.hotww and 1.0 < azz ^ 1-6. However when converted to signal 
strength it becomes fizz ~ Hww and 0.6 < Hzz,ww ^ 1-4. The best fit values for 
both cases are SM-like. 




0.8 1.0 1.2 1.4 1.6 ~ 1.8 '0.4 0.6 0.8 1.0 1.2 1.4 1.6 



azz fizz 

(a) (b) 

Figure 4: The 95%CL correlations between the of h — > WW* and h — > ZZ* decay rate 
ratios, sub-diagram (a), and signal strength, sub-diagram (b). Where the best fit point 
is marked by the dot/cross, and the blue (darker)/ brown(lighter) region is for using the 
LHC/SM data. 



5.3 Constraints on the gauge-Higgs Wilson coefficient by using the Higgs 
data plus EMD's and AMDM's 

Finally, we include 4 more data points, d e , d n , Aa e , and Aa^ into the global fitting. 
And the 95%CL results are shown in Fig.|5]. One can clearly see that when the 
electron and neutron EDM are included in the global fitting, the limits on the CP- 
odd Wilson coefficients are dramatically improved. Even giving extra, say ~ O(10), 
uncertainties to the 1-loop EDM estimations by using the effective operators, the 
inclusion of EDM's into the analysis still plays an important role to limit the CP 
properties of 125 GeV boson. More implications of including the EDM constraints 



will be discussed in Section 5.3 



Some interesting correlations emerge among a-yz,zz,ww, see Fig|| and Fig.^(a). 
Roughly, we observe that 

(a ww - 1) ~ 2(a Z z - 1) , 

{a lZ - 1) ~ 6.5 x 10 3 (a w - l) 2 ~ 2.57 x 10 4 (a zz - l) 2 . (5.6) 

These relations can be easily understood. Since the EDM's strongly constrain the 
CP-odd Wilson coefficients, the a's will be dominated by the CP-even Wilson coeffi- 
cients ci and C2. From Eq. ( |3.19| ), we see that a^zyvw.zz are more sensitive to c% than 



Ci, thus the approximate relations follow. However, when translated into the signal 
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Figure 5: The 95%CL allowed region in the Wilson coefficient space by using Higgs data, 
EDM, and AMDM. The blue(darker) and light brown (lighter) regions are for using the 
LHC data and SM predictions respectively. The best fit location is shown at the dot (cross) 
for LHC data (SM). The corresponding minimum x 2 = 15.89(9.63) for LHC data(SM). In 
sub diagram (b), only the SM-like C3 region is shown. 




1.04 



Figure 6: The 95%CL decay rate ratios correlation between: (a) h — > WW and h — > jZ 
, and (b) h — > ZZ and h — > 7Z. Where the best fit point is marked by the dot/cross, and 
the blue (darker)/ brown(lighter) region is for using the LHC/SM data. 



strength, only the relation between fizz and [iww is hold, see Fig.[7|(b). Therefore, 
we find that nzz ~ Hww and 0.6 < Hzz,ww ^ 1-4 is quite robust with or without 
taking EDM constraints into account. To test the relations given in Eq. (|5.6j ), better 
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than < 1% accuracy of determination the individual absolute Higgs decay width is 
required. This will be very challenging at LHC but could be done at the precision 
Higgs machines in the future. 



1.4 
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Figure 7: This plot is similar to Fig. 4 but with additional constraints from EDM's and 
AMDM's. 



6. Discussion and Conclusion 

6.1 Discriminate the degenerate solutions 

From the global fit of the latest LHC Higgs data, in the presence of effective gauge- 
Higgs operators, we obtain 7 99 ~ 1.0 ± 0.4. This introduces a degenerate solution to 
C3, 1 + 13.92C3 = ± a/1.0 ± 0.4. One solution is SM like with C3 centers around zero. 
Another possible solution, C3 ~ —0.14, corresponds to that the NP gluon-gluon fusion 
amplitude equals minus two times the SM amplitude. We would like to point out 
that this degeneracy can be easily lifted once LHC has enough Higgs pair production 
data in the high luminosity phase. The Higgs pair production by gluon-gluon fusion 
in SM, G^(p a )G u (pb) —> h(p c )h(pd), has been analyzed in [47] and summarized in the 
review paper [36]. In SM, the two gluons must carry the same color and this process 
receives contributions from the triangle diagram and box diagram. Both diagrams 
yield the spin-0 form factor while the box diagram contributes additional spin-2 form 
factor. The SM amplitude is 



M(G a G b ->• H c H d 



Gpct s s 



s — m 



2 FtA^ + (F B A 0flu + GbA 2/j , u ) 

h 



b i 



(6-1) 
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where s = (p a +Pb) 2 , and the spin-0/2 form factors Aq U (oc S^/A^ are given by 

(6.2) 

2 Oft • Pc)p2p£ _ 2(Pa-Pc)PbPc + 2 PcPc 
PriPa-Pb) PriPa-Pb) Pt 

For more details, see [47]. The operator 0% provides an additional contribution to 
the spin-0 amplitude, and the amplitude square becomes 

+ \G B \ 2 ^j ■ (6.3) 

If taking A = ITeV, the exotic solution that c 3 ~ -0.14 gives -f^§rc 3 ~ +1.34 
inside the spin-0 amplitude square. Because of the cancelation between the contri- 
butions from Ft and F R 2 , this exotic solution generates a sizable deviation from the 
SM prediction. Therefore, this two-fold degeneracy could be resolved by the future 
Higgs pair production data. 

6.2 A UV complete model 

Here we discuss a UV complete toy model to illustrate the validation of our effective 
operator analysis. Let's assume there is an additional color octet scalar S whose 
SM quantum number is (8,2,1/2) [33] on top of the SM. The most general scalar 
potential is 

V = - (H^Hi) 2 + 2m 2 s TrS^Si + X^H^iS^ S 3 + \ 2 H^H-TiS^ 

+ (X 3 H^ i H^ j TrS i S j + X^TtS^SjS, + X 5 H^TrS^ j SiSj + h.c.) 
+X 6 TrS jl S i S jj S :j + XjTrS^SjS^Si + X 8 TtS^S 1 TtS^ Sj 
+X 9 TrS fl S ] TrS^S i + X 10 Tr SiSjTrS^S^ + A u Tr S '<£ , , (6.4) 

where i,j are the SU(2) indices, and ins is the mass of the color octet scalar. The 
custodial symmetry requires the following relations to be held: 2A 3 = A 2 ,2A 6 = 
2A 7 = An, Ag = Aio for the real couplings [33], and A 4 = Ag for the complex ones [48]. 
And the SM quarks can now couple to both H and S by: 

-C Y = y^Q Ll Hu Rj +yf 3 Q Ll Hd R] + Y^Q u S A T A u Rj +YgQ u S A T A d R3 + hx. , (6.5) 

where T A is the SU(3) generator, and A is the color index. The Yukawa couplings 
y and Y are in general complex. In the Fig.|8| we show how we can generate these 

2 In the large quark mass limit F T = | + 0{s/m^), Fb = — § + 0(s/toq), and Gb — 0(s/mg) 
[47]. 
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Figure 8: Some typical loop diagrams which generate the effective gauge-Higgs operators 
in the color octet scalar model. Where S, the red(thick) dash lines, are the color octet 
scalars, H is SM Higgs doublet, V and V are the SM vector bosons. Note there are many 
other ways to attach the V and V if gauge symmetry allowed. Diagram (a) is the leading 
contribution to the CP-even gauge-Higgs operators. The CP-odd gauge-Higgs operators 
begin at two- loop level, diagrams (b). 



effective gauge-Higgs operators. One typical diagram to produce CP-even operator 
is Fig||(a), which is at one-loop level. Based on a naive dimensional analysis, the 
Wilson coefficients are given [33] 

c 3 _ 3 c 2 _ 3 ci 
A 2 ~ 2 A 2 ~ 2A 2 

A 2 ~~ 487r 2 m| ' 

And the CP violating operators are generated at two-loop levels, see Fig.|8|(b). A ball 
park estimation gives 

Cl,2,12,3 Im[A%] rfi 7] 

A 2 ^(16vr 2 ) 2 m 2 " { ' 

The strength of c's are roughly two orders smaller than the CP-even Wilson coeffi- 
cients in this model. 

Since this model is UV complete, we are able to discuss the fermion EDM and 
AMDM at above the electroweak scale without encountering any divergence. The 
charged lepton (g-2) starts at two-loop level, see Fig.|](a,b) 3 . Both Fig.[|(a) and 
Fig.^(b) can be related to the diagram shown in Fig^(a) by substituting VEV(s) for 
either one or two of external H legs. And this is exactly what we have performed in 
the gauge-Higgs operator analysis for Aa e and Aa M . 

For a most general Yukawa coupling, Eq. flfi.op , SM quarks receive nonzero EDM 
at 2-loop level, see Fig.|](c,d) 4 . Both Fig.||(c) and Fig.|9|(d) are independent of 

3 The resulting 2-loop (g-2) can be read and translated from [49] where the similar diagrams due 
to exotic scalars have been considered. 

4 The formulas in [50], where EDM's are generated via the similar diagrams in SUSY models, 
can be easily translated for use in this Octet model. 



2Ai + A 2 
647r 2 m| 



(6.6) 
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Figure 9: Some typical leading loop diagrams which generate quark or lepton AMDM, (a) 
and (b), and quark EDM and cEDM, (c) and (d), in the color octet scalar model. Where 
S, the red(thick) dash lines, are the color octet scalars, and / is the SM fermion, either 
lepton or quark. Note there are many other ways to attach the external photon if gauge 
symmetry allowed. 

Fig.|8|(b), and our gauge-Higgs operator estimation for quark EDM are subleading. 
On the other hand, there are two kinds contributions to SM lepton EDM: (1) by con- 
necting Fig.|8](b) to the lepton line by either SM Higgs or gauge bosons, these 3-loop 
diagrams have been taken care of by our gauge-Higgs operator analysis, and (2) by 
joining the two external quarks to form a loop in Fig.|9](c) and Fig.[5](d) and connect 
the resulting bulb to the lepton lines by either SM Higgs or gauge bosons. These are 
4-loop diagrams. So our gauge-Higgs operator estimation for lepton EDM are indeed 
the leading contribution. Nevertheless, in the case that Yukawa coupling are most 
general, one needs to incorporate the Fermion EDM/ AMDM effective operators in 
the analysis, see [51] for an earlier study on the interplay of the 4-fermi operators 
and the electric and magnetic diploe operators. 

However, to avoid the FCNC and other phenomenological problems, it is a com- 
mon practice to assume that Y^ D = (3 u / D y^ D , where (3 U / D are real. In that case, 
both Yukawa are real in the fermion mass basis, and the contribution from Fig.^|(c) 
and Fig.|9|(d) vanish. Then our EDM estimations based on the gauge-Higgs operators 
become the leading one. This octet model with MFV is a good example for our EDM 
and g-2 analysis to work based on the gauge-Higgs operators. 

Assuming that A = m s = ITeV, our numerical gives —0.75 < c x < 0.75, 
-1-0 < c 2 < 1.0, -1.0 < ci 2 < 1.0, -3 < ci < 4, -2.5 < c 2 < 3, -0.25 < c 12 < 0.25 
(from S parameter ), and —0.02 < C3 < 0.015, if we take the SM-like solution, then 

-118 < A 2 < 118, 
-12.6 <2Ai + A 2 < 9.5, 

I A 4 1 1 sin 6> 4 1 < 137, (6.8) 

where we write A4 = lA^e* 6 * 4 . In addition to the positivity conditions for the scalar 
potential, Eq.( |6.4j ), the above limits provide further nontrivial constrains on the 
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parameter space of this model. 



6.3 Predictions for the CP-odd decays 

The spin and parity of the 125 GeV boson has been largely determined to be J p = + 
from the analysis of the polar angular distribution in the h — >■ 2j mode [8, 16] (for 
spin) and the polarization correlation in the h — > ZZ — > 41 decay [9, 17] (for parity). 
However, the current spin-parity analysis are all based on limited number of events 
and simple working assumptions, namely the boson is assumed to be either purely 
scalar, pseudoscalar, or spin-2 object ect. But we have to keep in mind that even 
the 125 GeV boson is an elementary scalar, the CP violation decays can still be 
generated by quantum corrections. For a given model, whether we have the required 
experimental sensitivity to detect the CP-odd composition is the question. 

If at low energy the new physics can be indeed described by the gauge-Higgs 
operators alone, our numerical study predicts interesting model-independent rela- 
tions among the Higgs CP-odd decay modes. As has been discussed in Sec. 4, also 
from Fig.[5](e,f), the three CP-odd Wilson coefficients follow a linear relation on a 
line segment with end points at c\ 2 ~ 0.56/1.0 for using the SM predictions/ LHC 
data as input. For A = ITeV, we have (c.f. Eq.( p,19| )) 

< ~3.215 2 12 , 
a% -1.92 4, 
af z ~4.3x 10" 5 4, 

a^-l.lSxlO- 4 ^. (6.9) 

Note that both h — > WW and h — > ZZ channels are not sensitive to the presence of 
the CP-odd gauge-Higgs operators. 

We are interested in the CP violation fraction, CPVij = YqpjY % \ which is more 



relevant for the parity determination. In Fig.|l0|, we show the 95%CL CP violation 
fraction correlations between: (a) h — > ZZ and h — > 77, and (b) h — > WW and 
h — > 77. Note the possible large CPV fractions due to the gauge-Higgs effective 
operators in the h — > 77 and h — > ^Z modes, which are not constrained at all by the 
global fit. 

Apparently, it is not likely to probe the CP violation fractions, at the level of < 
( a few )xl0~ 5 , in the h — > WW, ZZ channels at LHC 5 . This agrees with the current 
parity determination results [9,17]. And hypothetically, any future observation of 
~ a few % level CP violation fraction in the h — > ZZ mode will definitely indicate 
the existence of new physics beyond the gauge-Higgs sector and additional effective 
operators must be included. 

For making use of h — > 27, 7Z channels to probe the CP-odd components of 
the 125 GeV boson, we need to study the spin correlations. This could be neatly 



3 See [52] for other proposals of measuring the CP nature of Higgs at LHC. 
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Figure 10: The 95%CL CP violation fraction correlations between: (a) h —> ZZ and 
h — > 77, and (b) h — > WW and h — > 77. The constraints are from EDM's, AMDM's, and 
LHC data (blue, darker) or SM predictions (brown, lighter). 



done in the future ILC and Higgs factory by the angular distributions of the final 
state fermions in h — > (7*7*), {l*Z) — > 4f. For that we have a model-independent 



prediction: qjp : 0^ ~ 1.67 : 1. However, as shown in Fig.|TT], when the above 
relation is converted into the CPV fraction, the relation is smeared but more or less 
follows a linear relation. From the plot, the CPV fraction correlation between these 
two modes is roughly 

CPV 1Z ~ (0.65 ± 0.35)CPV 77 (68%CL). (6.10) 

However, the slop becomes (0.92 ± 0.7) for the 95%CL fit. Although not being very 
predictive, this prediction can be checked once the experimental sensitivities meet 
the SM theoretical uncertainties and the CPV fractions in 77 and jZ modes are 
measured in the future. 

Nevertheless, from our numerical study, we found potentially large CP violating 
compositions in the h — > 77, 7Z decays. 

6.4 Summary 

We have studied the new physics associated with the 125 GeV boson where the NP 
is assumed to be characterized by the effective gauge-Higgs operators. By global 
fitting we studied the correlations among various Higgs to di-boson decay modes. In 
addition to the updated Higgs measurements accumulated at LHC, we also considered 
the case when the experimental sensitivities are compatible with the SM theoretic 
uncertainties. We found that there is plenty of room for new physics to hide in the 
shadow of the inherent SM theoretical uncertainties. From the global fit, we found 
a robust prediction that fizz — Hww and 0.6 < (J>zz,ww ^ 1-4 at 95%CL, see Figs|| 
and 0. This can be tested at LHC with more data in the near future. 
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Figure 11: (a) The 68%CL and (b) the 95%CL CP violation fraction correlations between 
h — > 77 and h — > jZ. The constraints are from EDM's, AMDM's, and the SM predictions. 

Moreover, we take into account the constraints of EDM's and AMDM's under 
the assumption that the gauge-Higgs operators give the dominate contributions to 
them. We found very strong constraints on the CP-odd Wilson coefficients and our 
numerical indicated that the CP violating hVV interactions in the h — > WW and 
h — > ZZ modes are too small, ~ O(10~ 5 ) at most, to be detected at LHC, see Fig.|T0[ 
However, the CP-odd fraction in the ^77 and hjZ channels could be as large as 0(1). 

Two more intriguing relations are predicted if the EDM constraints are included: 
(1) (otww — 1) ~ 2(azz — 1), and (2) (a 7 z — 1) ~ 6500(aww — I) 2 - These predictions 
could be tested at the future precision Higgs machines. 

The caveats: we should keep in mind that our results are valid only for the 
class of NP which manifests itself in the form of the gauge-Higgs operators discussed 
in this work. So our results are not applicable if (1) existence of any light exotic 
degree of freedom below the electroweak scale which interacts with the Higgs boson, 
or (2) the NP is beyond the realm which can be largely charted by the gauge-Higgs 
operators. On the other hand, any violation of our above predictions will indicate 
that the new physics must go beyond the gauge-Higgs sector. 
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A. Decay width of h ->• WW*, ZZ* 

Details of calculating the width of h — > VV* (V = W,Z) are collected in this section. 
We consider the case that the virtute gauge boson V* subsequently decays into two 
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fermions / and /' with their momentums labeled as: h(P) — > V(pi) + /(P2) + f'iPs)- 
In the center of mass frame of / and /', the three-body decay width is given by the 
following phase integral 

I r(M-mi) 2 rS2+ 



1 

s 2 ± = m\ + m 3 + — 



(s — Si — m 1 )(si + m 3 — m 2 ) ± A 2 (si,m 3 , m 2 )X 2 (si, M , m x ) 

A(a, 6, c)=a 2 + b 2 + c 2 - 2ab - 2ac - 2bc , (A.2) 

where the kinematical variables are defined as: s = P 2 , p A = (p 2 + Ps), s\ = p 2 , 
s 2 — (Pi + P3) 2 , S3 = (pi +P2) 2 , and Aii and Ai2 are the CP-even and CP-odd 
amplitudes respectively The phase space integral is largely simplified when the final 
state fermions are massless which is a good approximation in our case. 

Beside the SM contribution, the CP-even amplitude also receives the contribu- 
tions from Oi and 2 and 

l-Mil 2 = ^M±^l x \(m 2 v AV2G F - a 2 si) (2 Pl • p 2 pi ■ p 3 - p\p 2 ■ p 3 ) 
\ s i m v) 

+(m 2 v 2(V2G F )* - a 2 pi ■ p A ) 2 2p 2 • p 3 ] , (A.3) 

where a' w = 2c 2 vg 2 / 'A 2 , a' z = 2a A vg 2 / A 2 , A"/ the color factor of the final state 
fermions, and (fiv,fiu) = (g 2 /2c w )(I 3 -2Q f s 2 v , -J 3 ) for Z° and g v = -g A = g 2 /2\/2 
for W ± . And Ai 2 is completely beyond SM, 

\M 2 1 2 = 8(^ + ^) ( — J x [(pi -pa)(p3 -pm -P2+P2-P4P1 - Pa) 

\ s l m V / 

-p\p\ ■ P2P1 ■ Ps - P 2 iP2 ■ P4Ps ■ Pa] , (A.4) 

where d' w = 2c 2 vg 2 /A 2 , and d' z = 2a 4 vg 2 /A 2 . Summing up all massless final states, 
the corresponding CP-even decay widths are: 

Tf = 9 ^^ R 1 (m w /rn h ,a' w ) (for W) , 

rf = 3G \f Z 3 mh S z Ri(m z /m h , a' z ) (for Z) , (A.5) 

with 5^ = (7 — y s w + the phase space factor i?i will be given in below. 

The Sz factor is different from [36] and differs by a factor of 3 for W case, but our 
results are consistent with [53]. 

For the CP-odd part, the decay widths for W and Z bosons are 

r w ~2 9m h9%v 2 2 ^G 2 F m^ w m h ?>2m 2 h m 2 w 

V2 = ~327r^A^ = 87T^ A 1 ^ ' ' 

Tf = al5 z j^f ^^i? 2 (e) . (A.6) 
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The two functions R±, R 2 can be expressed in the form of integrals, 

(l-c) 2 ri2+ 1 



r(i-e)~ p 

Ii':(<-<i\ ) = e 2 / d5i / ... ^ . , , 

</0 •/S; 



(si - e ) 

2 ' 



R 2 (e) = I dst I dh—_^ 2 



S2- 



1(1 _ e 2 _ gl) 2 _ l £ 2 5i _ 1 ( ~ 2 _ £ 2 )(1 _ gi _ 



with the relevant parameters introduced as e = s± = s 2 = and s 2 ± = 

|(1 + e 2 — si) ± |a/[(1 + e) 2 — 5i][(l — e) 2 — §1]. The analytical form for Ri >2 can be 
obtained: 



Ri(e,a) = 



' „ (l-e 2 ) v / 4e 2 ^T _ ' 

A + 5 arctan h C In e 

1 — 3e 2 



(A.7) 



72e 2 

A = (1 - e 2 ) [36ae 2 (-5 + 9e 2 ) + 6(2 - 13e 2 + 47e 4 ) + a 2 (17 - 82e 2 + 89e 4 )] , 

B = J = [6(e 2 - 8e 4 + 20e 6 ) + 12a(e 2 - 8e 4 + 14e 6 ) + a 2 (-l + lie 2 - 40e 4 + 54e 6 )] , 

C = 6 [6(e 2 - 6e 4 + 4e 6 ) + 12a(e 2 - 6e 4 + 2e 6 ) + a 2 (-l + 9e 2 - 30e 4 + 6e 6 )] , 

R2(e) = { ^p- (-17 + 64e 2 + e 4 ) + ^ (1 + 6e 4 - 9e 2 ) 

(7e 2 -l) (e 2 - l)^ 2 - 1 lk x 

+ - V4e 2 -1 arctan ^ , (A.8) 

6 3e z — 1 



where 



^c 2 , for W 



- K 2 1 a 4 , for Z 



One can check that when a' v = 0, the SM result is recovered, and 6R1 = F(e) as 
defined in [53]. 
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